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The F-GAMMA program is among the most comprehensive programs that aim at understanding 
the physics in active galactic nuclei through the multi-frequency monitoring of Fermi blazars. Here 
we discuss monthly sampled broad-band radio spectra (2.6 - 142 GHz). Two different studies are 
presented, (a) We discuss that the variability patterns traced can be classified into two classes: (1) 
to those showing intense spectral-evolution and (2) those showing a self-similar quasi-achromatic 
behaviour. We show that a simple two-component model can very well reproduce the observed phe- 
nomenologies. (b) We present the cm-to-mm behaviour of three 7-ray bright Narrow Line Seyfert 
1 galaxies over time spans varying between ~ 1.5 and 3 years and compare their variability charac- 
teristics with typical blazars. 



I. BLAZAR VARIABILITY 

Among the most prominent characteristics of 
blazars, is their intense variability at all energy bands 
from radio to 7-rays and TeV energies. The mecha- 
nism producing variability itself depends on the en- 
ergy band it refers to and has been long debated. The 
"Shock-in- Jet" model suggested by [1], is broadly ac- 
cepted to explain the variability in terms of shocks 
propagating in the jet. The basic assumption is that 
changes in the particle injection rate, the magnetic 
field and the bulk Lorentz factor taking place at the 
jet base, cause the formation of shocks which conse- 
quently go through first Compton, then synchrotron 
and finally adiabatic energy loss phases. An alterna- 
tive model, the "Internal Shock Model" proposed by 
Q, suggests that energy is channeled into the jet in 
an intermittent way. "Plasma shells", may have dif- 
ferent bulk Lorentz factors and masses, so that faster 
shells can catch up with slower ones, collide with them 
and create relativistic shocks. The shock accelerates 
electrons to relativistic energies so that they can emit 
synchrotron and inverse Compton radiation. Other 
models explain the variability geometrically. Villata 
and Raiteri [3] for instance, suggested that orbital mo- 
tion and jet precession, possibly caused by a binary 
black hole system, can produce helical jet morpholo- 
gies which may then cause brightness and spectral 
variability through changes in the Doppler factor. 

Studies of the variability of the Spectral Energy 
Distribution (SED), ideally based on simultaneous 
datasets, probe the physics driving the energy pro- 



duction and dissipation in these systems (e.g. [4|, [5|) 
and can readily reveal some physical parameters such 
as the brightness temperatures, Doppler factors and 
the jet viewing angles. 

The current report deals with the variability charac- 
teristics of AGNs as they are observed in radio bands 
at cm to mm wavelengths and has a twofold character. 
First, we present the phenomenological classification 
of the variability patterns followed by selected Fermi 
blazars and propose a toy-model which can physically 
explain this classification. Secondly, we present the ra- 
dio behaviour of three Narrow Line Seyfert 1 galaxies 
(hereafter NLSyl) which have been detected in 7-rays 
by Fermi/LAT . 



II. THE F-GAMMA PROGRAM 

The potential of the LAT, on-board the Fermi satel- 
lite [6j, in the study of AGNs, could be attributed 
mostly to (a) the covered energy range (from 20 MeV 
to >300GeV) and (b) the pace at which the 47r 7- 
ray sky is being scanned. Given the broadness of the 
blazar spectral energy distribution (SED), the full po- 
tential of such densely sampled 7-ray light curves can 
be explored only when combined with multi-frequency 
data. Radio monitoring programs are essential in this 
task since they trace the behaviour of the relativistic 
jets where both primary and secondary emission pro- 
cesses take place (i.e. synchrotron and inverse Comp- 
ton processes). 

One of the most comprehensive radio monitoring 
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efforts, the F- GAMMA program, [7-9] was initiated 
in 2007 utilising the Effelsberg 100-m, the IRAM 30- 
m and the APEX 12-m telescopes for the monthly 
monitoring of ~ 60 Fermi blazars at 12 frequencies 
between 2.6 and 345 GHz. The millimetre observa- 
tions are closely coordinated with the more general 
flux monitoring conducted by IRAM at the 30-m tele- 
scope. The Effelsberg 100-m telescope is equipped 
with circularly polarised feeds, while the IRAM 30- 
m with linearly polarised ones. Details are given 
elsewhere (Fuhrmann et al. in prep., Angelakis et 
al. in prep., Nestoras et al. in prep.). Measure- 
ments at 4.85 GHz, 10.45 GHz, 32.0 GHz, 86.24 GHz 
and 142 GHz are done differentially either by using 
multi-feed systems, or, at IRAM 30-m, by wobbler 
switching. On average, the time needed for observing 
an entire spectrum of any given source at Effelsberg 
alone is of the order of 35 minutes, while at IRAM, 
roughly 2 minutes. The combined spectra (Effelsberg 
and IRAM) are observed quasi-simultaneously within 



approximately one week. That is, neither the single- 
facility spectra nor the combined ones are likely to be 
affected by source variability. In the current study 
only data collected until June 2011, have been used. 



III. BROAD-BAND RADIO SPECTRA 
VARIABILITY 



The data product of the F-GAMMA program is 
monthly sampled broad-band radio spectra covering 
a frequency range of 2 orders of magnitude. In the 
following we discuss studies based only on Effelsberg 
and IRAM data i.e. 2.6 - 142 GHz and present a toy- 
model for the interpretation of the small number of 
phenomenological types of variability patterns. 
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FIG. 10: The assumed two- component system. The dif- 
ferent variability types can be reproduced with the ap- 
propriate modulation of the relative position and relative 
broadness of the band-pass denoted by the grey shaded 
areas. 



A. Phenomenological Classification 

Angelakis et al. PH [TTJ] showed that the variabil- 
ity patterns that the source spectra follow as a func- 
tion of time, fall in a small number of phenomeno- 
logical classes the prototypes of which are shown in 
Figures [T] - [9] A detailed description of the charac- 
teristics of each type are given in [10|, What is 
readily obvious in this analysis and can also be seen 
in Figures [T] - [9j is that the patterns can be basically 
classified in (a) those characterised by spectral evolu- 
tion (Figures [T] - [7j) and (b) those that show a convex 
spectrum which varies self-similarly following a close- 
to- achromatic evolutionary path (Figures [8] and [9j). 



B. Interpretation of the Observed 
Phenomenologies 

As it is argued in [10], the spectral-evolution dom- 
inated types can naturally be explained by a simple 
two-component system, as it is observed by other pro- 
grams as well [e.g.[l2|. It is shown that a system com- 
posed of a steep spectrum component attributed to a 
large scale jet populated with a varying spectral com- 
ponent, can easily reproduce the observed patterns 
(Figure [10]). The spectral component is assumed to be 
following the [l| model. Calculations showed that with 
only a minor coverage of the parameter space that de- 
termines the characteristics of those two elements (e.g. 
intrinsic source luminosity, flare dominance, turnover 
frequency of quiescent spectrum), is enough to repro- 
duce the observed variability. This discussion refers 
only to these types that are characterised by spectral 



evolution. The "achromatic" behaviour must be at- 
tributed to different mechanism as it will be discussed 
elsewhere. 



C. Variability Properties 

The different classes discussed earlier have some 
very fundamental differences. The most obvious one 
being the increase in the prominence of a quiescent 
large scale jet as compared to the flaring event with 
increasing type. That should immediately imprint dif- 
ferences in the distribution of the variability parame- 
ters for types that have different quiescent spectrum 
prominence. 

In order to investigate that, we have performed 
standard variability analysis. For each light curve 
characteristic variability time scales have been esti- 
mated using a structure function analysis described 
by [l3[. These time scales are subsequently used to 
compute the involved emitting region sizes (assuming 
the light- travel-time argument) and hence the vari- 
ability brightness temperatures, as: 

r b[ i<i = «.io.»A5( z7 ^ ; ) 2 (1) 

where AS: the flux density variation in Jy 
A: the observing wavelength in cm 
D^: the luminosity distance in Mpc 
At: the characteristic time scale in days 

Assuming a typical value of 5 • 10 10 K for the up- 
per limit in the brightness tem per ature on the basis 
of equipartition arguments [14], |l5| , we can compute 
also the variability equipartition Doppler factor, as: 

L> var = (1 + z) • (T b /5 • 10 10 K)^, where a is the 
spectral index (S oc v a ). A typical value of a — —0.7 
has been used. 

In Figure [11] we show the variability brightness tem- 
perature (Eq. [T]) distribution at three characteristic 
frequencies, 4.85, 14.6 and 32 GHz, for three groups 
of sources separately. One group is made of sources of 
type 3-4 which show the clear presence of a quiescent 
spectrum (qs); sources of type 1-2 comprise the second 
group and show no evidence for the presence of such a 
spectrum. The last group is made of sources that be- 
long to type 5. In the lower row of Figure [11] we show 
the same distributions for the variability equipartition 
Doppler factors, D. 

From these plots it can be seen that at all bands 
there is a separation between the distribution of the 
first group when compared with the other two ones. 
This could be explained as follows; for sources of type 
1-2 or 5 the dominance of flares relative to the qui- 
escent spectrum is significant. Given the fact that 
the brightness temperature and the corresponding 
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FIG. 11: Variability brightness temperatures and the inferred equipartition Doppler factor distributions. Red denotes 
sources of variability type 1-2 that show no clue of quiescent spectrum. Grey denotes sources where there is an evident 
quiescent spectrum and blue are sources of type 5. Left: 4.85 GHz, middle: 14.6 GHz, right: 32 GHz. 



Doppler factor are inferred from the variability prop- 
erties, it is naturally expected that sources where the 
quiescent spectrum is absent, show more pronounced 
variability characteristics. At the highest frequencies 
this separation must disappear due to the absence of 
quiescent spectrum. 



IV. JET EMISSION FROM NLSY1S 

Narrow Line Seyfert 1 galaxies are an AGN sub- 
class that have been argued to show comparably low 
masses and high accretion rates [e.g. [l6l4l9l|. O nly a 
mere 7 % of them appears to be radio loud |20j . 

The early Fermi discovery of 7 -ray emission from 
a small number of NLSyls [2l|, |22| ) came as sur- 
prise; until then the only 7-ray bright classes were 
thought to be blazars (i.e. FSRQs and BL Lac ob- 
jects) and radio galaxies (for a review, see [23]). The 
discovery of 7-ray emission from NLSyls not only in- 
troduced a new class of 7-emitting AGNs but also 
revolutionised the well spread belief that jets are as- 
sociated (chiefly) with large elliptical galaxies. The 
multi-wavelength campaign between March and July 
2009 following the discovery of 7-ray emission from 
PMN 0948+0022 showed that the source was exhibit- 
ing a spectral behaviour typical of a relativist ic jet 

Responding to the Fermi detection the F-GAMMA 
team initiated a dedicated program, to understand 
their cm to mm behaviour [26|. The three sources 



TABLE I: The NLSyl monitored by the F-GAMMA pro- 



gram. 



Source 


z 


Art 


/* 


(S4.85) 








(months) 


(Jy) 


1H 0323+342 


0.061 


18 


1.0 


0.40 


PMN J0948+0022 


0.585 


42 


0.9 


0.23 


PKS 1502+036 


0.409 


18 


0.9 


0.52 



^ N: number of observed epochs until February 2012 



/: mean sampling 



that are monitored are listed in Table HI In the follow- 
ing we discuss first results of this monitoring. 

A. Radio Variability 

In Figure [12] we show the variability light curves as 
well as the broad-band spectra observed within the 
context of the F-GAMMA program. For 1H 0323+342 
and PMN J0948+0022 light curves between 2.6 and 
142 GHz have been constructed while the steepen- 
ing of the spectrum above 32 GHz in the case of 
PKS 1502+036, does not allow IRAM observations at 
short-mm bands. 

1H 0323+342: As it can be seen in Figure H2J 
its radio spectrum displays a quiescent part reaching 
up to roughly 10 GHz beyond which a high frequency 
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FIG. 12: Monthly sampled spectra and light curves for the three monitored NLSyls. Upper row: the light curves at all 
F-GAMMA frequencies between 2.64 and 142 GHz. Lower row: the broad-band F-GAMMA spectra. 



component (hereafter HFC) appears. Its mean 
spectral index is of the order of —0.23 (calculated 
between 2.6 and 10 GHz), not very different from 
a typical value of —0.5 (S ex v a ). The HFC shows 
intense spectral evolution which gradually shifts 
its peak progressively towards the steep spectrum 
component. The pace of the evolution is remarkably 
fast causing significant displacements within the 
month that typically separates two observations. 
The variability pattern according to the classification 
discussed in Sect. IIII Al is of type 4. Its light curves 
although only shortly sampled, shows a collection 
of events more prominent at higher frequencies 
and with cross-band lags indicative of the spectral 
evolution superimposed on a long term decreasing 
trend. The modulation index m (m[%\ = 100 x ) 
at 4.85 GHz is, 777,4.85 ~ 10%, while at 14.6 and 
32 GHz it is, 17% and 27%, respectively. The 
Structure Function analysis applied on the source 
light curves revealed characteristic times scales of 
the order of 60 days which implies a variability 
brightness temperature of 10 12 K at 4.85 GHz and 
2 x 10 11 K at 14.6 GHz. The corresponding equipar- 
tition Doppler factors are 2.4 and 1.5 respectively, 
placing the source in the lower part of the Doppler 
factor distribution shown in Figure [12] Interest- 
ingly, the Doppler factor calculated from fitting the 
Spectral Energy Distribution (SED) [2l[ is around 17. 

PMN J0948+0022: It is evident from the light 
curves and the spectra shown in Figure [12] that 
the available dataset for PMN J0948+0022 is much 



richer than that for the other two NLSyls. The 
spectrum appears mostly inverted representative of 
variability type 1. The spectral index below 10 GHz 
ranges between marginally steep or flat (« —0.1) to 
highly inverted reaching values of +1.0. Its evolution 
is exceptionally dynamic with significant evolution 
happening even within one month. It appears that 
for this source a sampling of two weeks would be 
necessary. Its light curve shows at least 4 prominent 
events which emerge with time lags at different 
bands. At the lowest frequencies the events are 
barely seen. Yet, the modulation index shows a 
monotonic increase with frequency. At 2.6, 4.85, 
14.6, 32 and 142 GHz the modulation index is 22, 
29, 35, 37 and 38%, respectively. The standard 
variability analysis reveals brightness temperatures 
of 8 x 10 12 K at 4.85 GHz, 2 x 10 12 K at 14.6 GHz 
and 1.5 x 10 n K at 32 GHz which would imply rather 
typical Doppler factors of 6, 4 and 2 respectively, as 
it can be seen from the Doppler factor distribution 
plots in Figure [11] The SED modelling gives Doppler 
factors that vary between 10 and 20 with the latter 
being observed during the outburst of July 2010 [27j. 

PKS 1502+036: This source shows a variabil- 
ity behaviour similar to type lb although more 
epochs are needed for a definite classification. Its 
spectrum is highly variable displaying periods of 
convex shape. Its low-band part (y < 8 GHz) varies 
between flat and highly inverted, while the higher 
frequencies (10 GHz < v) can show spectral index 
as steep as —0.5. The light curve shows at least 2 



eConf C1111101 



6 



2011 Fermi & Jansky: Our Evolving Understanding of AGN, St Michaels, MD, Nov. 10-12 



events better seen at higher frequencies. The high 
frequency cut-off of the spectrum prohibits IRAM 
monitoring. The typical time scales identified here 
are of the order of 60 - 80 days. At 4.85 GHz the 
brightness temperature is 2 x 10 13 K and at 14.6 GHz 
it is 3 x 10 12 K implying Doppler factors of 7 and 4, 
respectively while the SED fitting gives a value of 18. 



modelling. In any case, it seems premature to draw 
any final conclusions for the behaviour of this limited 
sample of NLSyls. Longer time baselines will allow 
more accurate estimates of the brightness tempera- 
tures and all the variability parameters. 



B. Conclusions 

From the above discussion it becomes obvious that 
the three monitored NLSyls show a very typical 
blazar-like behaviour. That is, highly variable spec- 
tra caused by the presence of prominent evolving high 
frequency spectral components. The variability hap- 
pens at interestingly fast pace with the mean number 
of events per unit time being clearly larger than that 
of the rest of the F-GAMMA targets. As an exam- 
ple, PMN J0948+0022 shows, on average, 1.4 flaring 
events per year as compared to less than ~ 1 event per 
year for other typical F-GAMMA blazars. The vari- 
ability brightness temperature is relatively high with 
respect to the distribution of the whole sample in the 
cases of PMN J0948+0022 and PKS 1502+036. The 
opposite is the case for 1H 0323+342. The derived 
Doppler factors are lower limits and are systemati- 
cally lower than the values obtained from the SED 
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